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We study interaction effects on transport through a small metallic cluster connected to two fer- 
romagnetic leads (a single-electron spin- valve transistor) in the "orthodox model" for the Coulomb 
blockade. The non-local exchange between the spin accumulation on the island and the ferromag- 
netic leads is shown to affect the transport properties such as the electric current and spin-transfer 
torque as a function of the magnetic configuration, gate voltage, and applied magnetic field. 



Magnetoelectronics is a contender to fulfill the tech- 
nological need for faster and smaller memory and sens- 
ing devices. The drive into the nanometer regime brings 
about an increasing importance of electron-electron inter- 
action effects. Small metallic clusters (islands) that are 
electrically contacted to metallic leads by tunnel junc- 
tions and capacitively coupled to a gate electrode can 
behave as "single-electron transistors (SETs)". In the 
Coulomb-blockade regime, the charging energy needed to 
change the electron number on the island by one exceeds 
the thermal energy and transport can be controlled on 
the level of the elementary charged In a spin-valve SET 
(SV-SET), the contacts to the cluster consist of ferromag- 
netic metals (F). We focus here on F|N|F structures with 
normal metal (N) islands (see Fig. da))^ since these 
"spin valves" display giant magnetoresistance and spin- 
current induced magnetization reversal. 3 Other combina- 
tions such as F|F|F— , N|N|F£ or F|F|N&£ are of interest 
as well. 

Several theoretical studies have been devoted to the 
binary magnetoresistance (MR) of SV-SETs, i.e. the 
difference in the electric resistance between parallel 
and antiparallel configurations of the magnetization 
directions&S Interaction effects in magnetic devices have 
been studied as well for spin valves with a Luttinger liq- 
uid island 1 * 3 and for single-level quantum dotsii*i2iH with 
non-collinear magnetic configurations. 

A necessary condition for a significant MR in F|N|F 
structures is a spin accumulation on the normal metal 
island, viz. a sufficiently long spin-flip relaxation time 
Tgf. Seneor et alH have measured the MR of SV-SETs 
with gold islands with a r s f ~ 800 ps. A r S £ in the mi- 
crosecond regime has been reported for Co nanoclustersA 
The long spin-flip times in small clusters is not yet fully 
understood; it might simply be due to the probability of 
finding zero impurities in a given small cluster. 

In this Rapid Communication we elucidate the new as- 
pects of electron transport related to non-collinear mag- 
netization directions in metallic SV-SETs. It turns out 
that an effective exchange effect between the spin accu- 
mulation and the magnetizations has to be taken into 
account. 

We take the junction resistances sufficiently larger than 
the quantum resistance Rq = h/e 2 so that the Coulomb 
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FIG. 1: The spin-valve single-electron transistor, (b) The 
tunneling rates between the leads and the cluster depend on 
the spin accumulation s. 



blockade can be treated by lowest-order perturbation the- 
ory. We furthermore disregard the size quantization of 
states in the clusters, thus adopting the well-established 
"orthodox" modeli In our model system (cf. Fig. |IJb)), 
the ferromagnetic leads are treated as reservoirs with 
single-domain magnetization directions rh\ and fh%. Dis- 
regarding magnetic anisotropics, the relevant parameter 
is the angle 9 between the magnetizations. The capac- 
itances of the junctions are C\ and C%. The cluster is 
capacitively coupled to the gate, with capacitance Cq <C 

We assume a separation of time scales between the en- 
ergy relaxation that rapidly thermalizes injected charges 
and the slow spin relaxation. In this regime, the 
quasi-equilibrium excess spin s on the normal metal 
island is well-defined. In second quantization s — 

(h/2) J2kss> ( c L^ as/Cfc »') , where a is the vector of Pauli 
spin matrices and k, s denote the orbital and spin in- 
dices of the island states, respectively. This corresponds 
to a chemical potential difference (spin accumulation) 
A/i = 2S\s\/H, where S is the average single-particle 
energy separation (in terms of the static susceptibil- 



ity Xs '■ A/i = 2fxf 3 \s\ / [xsfi>))- Spin-flip relaxation is 
parametrized by the spin-flip time r s / or spin-flip conduc- 
tance G s f = e 2 1 (2r s /(5). For metals, 5 is much smaller 
than the thermal energy except for very small particles 
(diameter <5 nm). We restrict our attention here to a 
regime in which the bias energy is small compared to the 
thermal and the charging energies, but large compared to 
<5, so that the transport properties do not depend on the 
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energy-relaxation rate^ The island state is then charac- 
terized by the excess number of electrons n and net spin 
angular momentum s. A state distribution on the island 
governed by free-energy minimization under these spin 
and charge constraints can be used. 

The master equation for electron transport in the or- 
thodox model is determined by a tunneling Hamiltonian 
that can be treated by perturbation theory. We derive 
the appropriate Hamiltonian by collecting the leading 
terms in the transition probabilities. A crucial param- 
eter is the mixing conductance for a N|F junction 15 
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where n and m denote the transport channels in the 
normal metal, and r™ m and r" m are the correspond- 
ing spin-dependent reflection coefficients. The real 
part of the mixing conductance determines the spin- 
transfer torque that strives to align the magnetiza- 
tion to the spin accumulation^ In the limit of a tun- 
neling contact, ReG^ = (Gf + Gj_) /2, where G s = 

( e2 /h) J2 nm 0>mn — It's" 1 ! 2 ) i s the conventional tunnel- 
ing conductance for spin s. The torque is transferred by 
the electrons that tunnel through the contact and is in- 
cluded in the standard spin-dependent tunneling Hamil- 
tonian. The imaginary part of the mixing conductance 
can be interpreted as an effective field parallel to the 
magnetization direction of the ferromagnet. For a normal 
metal separated from a Stoner ferromagnet by a specular, 
rectangular barrier we can directly solve the Schrodinger 
equation. With Fermi momenta of the ferromagnet of 
/cft = 1.09A -1 , k F i = 0.42A" 1 (characteristic for Fe), 
a normal metal Fermi energy of e_F = 2.6 eV, a barrier 
height of 3 eV, and free electron masses (cf. Ref. 0) we 
find ImG Ti /(G T + G X ) = -0.26. This illustrates that, 
in contrast to metallic interfaces^ the imaginary part 
of the mixing conductance can be significant for tunnel 
junctions. This result is not sensitive to the width of 
the barrier but may depend strongly on material com- 
bination and interface morphology. For tunneling barri- 
ers made from magnetic insulators, 1 ' 7 Re G^ and ImG^ 
may become large compared to (Gj- + G{) , which results 
in different physics^ 

The Hamiltonian for the SV-SET reads: 



H = #* 
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where Hjq describes the electrons on the normal metal 
and includes the electrostatic interaction energy^ 
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For the two ferromagnetic leads (a = 1,2) Hp a — 
J2ks £ aks<^ a ks a aks and n denotes the excess electron num- 
ber on the island. The tunneling Hamiltonian for each 



contact reads H Ta = z~2 k qs T akqsa aks c aqs + h.c. Finally, 
the imaginary part of the mixing conductance gives rise 
to an effective exchange effect: 
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where A CXQ = -MlmG^/ (2e 2 ) . We note that the ef- 
fect of Eq. on the spin accumulation is identical to 
that of an external magnetic field applied in the direc- 
tion rh a . Such an effective exchange Hamiltonian has 
been introduced before for a Luttinger liquid attached to 
ferromagnetic leads piS Physically, electrons in N feel the 
ferromagnet through their tunneling tails that cause a 
spin-dependence of the reflection coefficients. The small 
spin splitting due to A ex in the ground state does not in- 
fluence the transport in the leading order of perturbation 
theory. 

Recently, the angular dependence of transport through 
spin valves has been studied for single-level quantum dot 
islands In these systems an effective field was found 
to act on unpaired quantum dot electron spins that is 
caused by virtual particle exchange with the leads in the 
Coulomb blockade. These correlations cause effects sim- 
ilar to those discussed here, but their physical origin is 
completely different. Eq. Q is caused by electron ex- 
change on a very fast time scale corresponding to the 
reciprocal Fermi energy and reflects the electronic struc- 
tures of the junctions, independent of the applied volt- 
ages and charging energies. In contrast, the correlation- 
mediated exchange is induced on time scales of the recip- 
rocal charging energy, changes sign with gate voltage and 
does not vanish for normal metal contacts*^ We find that 
also for classical islands the correlation exchange field can 
be of the same order as A cx , but since both effects can 
at least in principle be distinguished experimentally by 
gate voltage and temperature dependence, a more de- 
tailed discussion is deferred to a future publication. 

We introduce spin-dependent conductances for both 
junctions, G as = ire 2 p N p Fas T as /h (a = 1,2, s =1,1). 
Pn is the density of states at the Fermi level in the 
normal metal, and pFas the spin-dependent density of 
states in ferromagnet a. T s is proportional to the average 
tunneling probability over all channels for spin s, T s = 



\T„ 



The conductances are assumed to be con- 



stant within the energy interval of the charging energy, 
which is a safe assumption for metals. We introduce the 
total conductances G a = G a -\ + G Q |, the polarizations 

The tunneling rate for adding an electron through con- 
tact a (=1,2 = +, — ) in the considered regime where 
eV <C UbT then reads: 



nfl — *n+l 



(V,q,s)= C ^F(-E a (V,q)) 



(-E a (V,q))^(rn a .s), (5) 

where E a (V, q) = an a eV — e (q — e/2) / (Gi + G 2 ) is 
the electrostatic energy difference associated with the 
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tunneling of one electron into the cluster from lead a 
with q = —ne + CqVq the charge on the island and 
K a = (1/Ci + I/C2) 1 jC a . The other rates can be found 
analogously. We can also find the net spin current into 
the cluster due to tunneling events in which an electron 
is added to the island through contact a: 
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The dynamics of the spin accumulation is also affected by 
spin-flip scattering, the exchange effect discussed above 
and an external magnetic field B : 
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The master equation for the SV-SET is dp n /dt — 

- Pn (r»-™+ 1 + r™^"- 1 ) + Pn+1 r i+1 ^ n + Pn - 1 r n - 1 ^ n 

combined with ds/dt described above. Here, p n is the 
probability distribution for the number of electrons on 
the island and y n ^ n + 1 denotes the rate for adding an- 
other electron when the cluster has n excess electrons. 
We now focus on a quasi-stationary state dp n /dt = 0. 
When eV <C fc^T <C e 2 /2C the number of electrons 
on the island fluctuates between two values denoted by 
n = "0" and "1" . We can use detailed balance symmetry 
to find the stationary state, whence poT ^ 1 = piT 1 ^ . 
We find the conductance G: 



G(V G ,s)=G KS (V G ) 
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where G KS (V G ) = /?A mln GiG 2 / [2 (G x + G 2 ) sinh/?A min ] 
with A m j n = e {CqVg — e/2) / (d + C2)) describes one 
Coulomb blockade oscillation^ 

The spin accumulation is obtained from the condition 
ds/dt = 0. For a symmetric setup, in which the conduc- 
tance parameters that characterize the two contacts are 
equal, we find: 
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When the imaginary part of the mixing conductance van- 
ishes, the angular dependence of the conductance is pro- 
portional to a cosine. The result of Brataas et al. for 
a F|N|F spin valve without interaction^ is recovered by 
substituting /3A min / (2sinh/3A min ) by 1 (cf. Ref.0). 

In Fig. Eta) an example of the dependence of the con- 
ductance on the angle 9 is shown. When the magne- 
tizations are non-collinear, the exchange effect reduces 
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FIG. 2: Results for a SV-SET with G 2 = 2Gi/3, Pi = 0.8, 
P 2 = 0.7, Gsf = Gi/5, ImG^ = Gi/4, ImG^ = G 2 /4 and 
A m i n = 0. (a) Conductance as a function of the angle be- 
tween the magnetizations of the leads, (b) Spin accumulation 
components in the direction of mi (dotted), z xtni (dashed), 
and z (solid), (c) Spin transfer torque on ferromagnet f in 
the direction zx fh\ (dashed) and z (solid), (d) Conductance 
vs. magnetic field applied in the direction z Xmi for 8 — 
(dotted), 9 = ty/2 (dashed) and 6 = vr (solid). 



the spin accumulation and increases the conductance. In 
Fig. Htb), the spin accumulation as a function of angle 
illustrates that for non-collinear magnetizations, the non- 
local exchange pulls the spin accumulation vector out of 
the plane of the magnetizations. 

The spin current between the ferromagnets and the 
normal metal gives rise to a spin-transfer torque T a = 
—m a x (I s x fh a ) on ferromagnet a (see Fig. Etc)). I s is 
the net spin flow out of the ferromagnet, and is strongly 
modulated by the gate voltage. When the Coulomb 
interaction suppresses the current, the exchange effect 
becomes relatively more important. Fig. Hjd)) shows 
that, for non-collinear configurations, the exchange effect 
causes an asymmetry in the Hanle effect with respect to 
the sign of an applied external magnetic field. 

The curves in Fig. show the gate- voltage modulated 
conductance for a symmetric spin valve with half- metallic 
ferromagnetic leads (P = 1). Spin-flip is disregarded. 
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FIG. 3: Conductance (G) as a function of gate voltage for 
9 = (dotted), 9 = 0.6vr (dashed), and 6 = 0.9tt (solid). The 
ferromagnets are half-metallic and the conductance of both 
contacts is Gi. ImG\ l = ImG^ = Gi/4, G s} = 0. 

When the magnetizations are parallel (9 = 0), no spin 
accumulates on the island and the Coulomb-blockade os- 
cillation equals that of all-normal metal systems. When 
the angle 9 is increased, the conductance is suppressed by 



a counteracting spin accumulation. The exchange effect 
acts to reduce this spin accumulation for non-collincar 
configurations and can cause a local conductance mini- 
mum at A m ; n = (e.g. for 9 = 0.9 tt). Deep in the 
Coulomb blockade, a significant spin accumulation is pre- 
vented from building up and all curves converge. 

Interestingly the angular magnetoresistance for Lut- 
tinger liquids with ferromagnetic contacts 10 : looks very 
similar. In order to find experimental evidence for spin- 
charge separation it is therefore necessary to avoid spu- 
rious effects caused by the Coulomb blockade. 

In conclusion, we studied the transport characteristics 
for non-collincar spin valves in the Coulomb blockade 
regime. A non-local exchange interaction between the 
spin accumulation and the ferromagnets affects the con- 
ductance and spin-transfer torque as a function of the 
gate voltage. This might provide new possibilities to 
control charge and spin transport in nanoscale magne- 
toelectronic devices. 
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